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Abstract. Aims. We present a detailed study of the infrared line emission (1-200/^m) in the Herbig-Haro ob- 
ject HH54. Our database comprises: high- (R~ 9000) and low- (R~600) resolution spectroscopic data in the 
near-infrared band (1-2.5/nn); mid-infrared spectrophotometric images (5-12^im); and, far-IR (45-200^im, R~200) 
spectra acquired with the ISO satellite. As a result, we provide the detection of and the absolute fluxes for more 
than 60 molecular features (mainly from H2 in the near- and mid-infrared and from H2O, CO and OH in the 
far-infrared) and 23 ionic lines. 

Methods. The H2 lines, coming from levels from v=0 to v—4 have been interpreted in the context of a state-of-art 
shock code, whose output parameters are adopted as input to a Large Velocity Gradient computation in order to 
interpret the FIR emission of CO, H2O and OH. 

Results. The H2 emission can be interpreted as originating in either steady-state J-type shocks or in quasi-steady 
J-type shocks with magnetic precursor. However, our multi-species analysis shows that only a model of a J-type 
shock with magnetic precursor (v s h oc k = 18 kms -1 , iih=10 4 cm~ 3 , B=100/iG, age=400 yr) can account for both 
the observed H2 emission and the CO and H2O lines. Such a model predicts a H2O abundance of ~ 7 10~ 5 , in 
agreement with estimations from other shock models for outflows associated with low mass protostars. We can 
exclude the possibility that the observed atomic lines arise in the same shock as the molecular lines, and give 
arguments in favour of the presence of a further high-velocity, fully dissociative shock component in the region. 
Finally, in view of the forthcoming spectroscopic facilities on board of the Herschel satellite, we provide predictions 
for H2O lines considered to be the most suitable for diagnostic purposes. 

Key words, stars: circumstellar matter - Infrared: ISM - ISM: Herbig-Haro objects - ISM: individual objects: 
HH54 - ISM: jets and outflows 



1. Introduction 

It is now commonly accepted that the earliest stages of 
star formation are accompanied by powerful bipolar jets, 
which act to remove excess angular momentum from in- 
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falling gas (e.g. Bacliiller, 1996). The jets are generally 
highly collimated and often terminate in a curve-shaped 
structure (bow-shock); however, knots of shocked emis- 
sion, which may display a bow morphology, are also ob- 
served along the jet (e.g. Reipurth & Bally, 2001). The 
thermal energy associated with a shock is radiated away 
through the emission of lines from atomic, ionic and 
molecular species. Gas at temperatures of thousands of 
degrees Kelvin cools principally through ro-vibrational 
lines of molecular hydrogen, in the near-infrared, but also 
through forbidden transitions of abundant atomic and 
ionic species, such as [Fe 11] , [S 11] , [C 1] and [Ni]. Colder 
gas components, at hundreds of degrees Kelvin, cool via 
mid- and far-infrared molecular lines, particularly rota- 
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tional transitions of H 2 (at A < 28 /zm) and lines of other 
molecular species, such as H 2 0, CO and OH (in the far- 
infrared) . The relative contribution to the gas cooling de- 
pends strongly upon the physical conditions within the 
shock. In particular, the shock velocity and the conditions 
in the pre-shock gas (the degree of ionization, the strength 
of the magnetic field and the density) determine the struc- 
ture of the shock wave (Draine, 1980), which is usually 
classified as 'jump' (J-type), 'J-type with magnetic precur- 
sor' or 'continuous' (C-type) (e.g., Hollenbach & McKee 
1989). In general, J-type shocks produce a narrower profile 
and a higher peak temperature than C-type shocks while 
J-type shocks with magnetic precursor can be considered 
to be intermediate. J-type shocks can heat the gas up to 
temperatures of 10 5 K and can travel at up to hundreds 
of kms" 1 , while C-type shock velocities do not exceed w 
80 km s" 1 (Lc Bourlot et al. 2002, Flower et al., 2003). 
Classically, calculations of shock structure have been per- 
formed under the assumption of steady state. However, 
if this assumption is relaxed it has been shown that, in 
the presence of a sufficiently strong magnetic field, an ini- 
tially J-type discontinuity may evolve in time into a C- 
type shock wave (Smith & MacLow, 1997, Chieze, Pineau 
des Forets & Flower, 1998). The intermediary stage of the 
evolution is described by a (quasi-steady) J-type shock 
wave with magnetic precursor and can be considered as a 
J-type discontinuity introduced into a C-type shock pro- 
file at the point in the ion fluid flow time, U = J 1/vidz, 
that may be identified with the age of the shock (Lesaffre 
et al., 2004a). 

Although 2D and 3D hydrodynamical and MHD 
simulations of jet-driven molecular shocks have pro- 
gressed enormously over the last years, including some 
time-dependent H 2 chemistry, H ionization, and cooling 
(Downcs & Cabrit, 2003; Smith & Rosen, 2003), the effect 
of ion-neutral decoupling in magnetic precursors and the 
resulting complex shock chemistry have so far only been 
simulated in ID (Chieze et al., 1998; Lesaffre et al., 2004b). 
Lesaffre et al. (2004a) showed that truncated steady-state 
multi-fluid models provide a useful and efficient approx- 
imation to time-dependent simulations, while allowing a 
more detailed treatment of the chemistry and hence of the 
degree of ionization. We shall thus adopt this approach 
here. 

Since different shock types give rise to very different ex- 
citation conditions, the interpretation of multiple species 
observed over a wide range in wavelength provides strong 
constraints for a shock model. Multi-wavelength studies 
have already proved useful: for example, Smith, Froebrich 
& Eisloffcl (2003) employed an analysis of the near- to far- 
infrared molecular emission observed towards Cepheus E 
in order to determine both the shock type and geome- 
try. The main limitation of the multi- wavelength/multi- 
species approach arises from the poor sensitivity and, 
more importantly, the poor spatial resolution achieved in 
the mid- and far-infrared by space-borne instrumentation, 
which often prevents the resolution of structures along jets 
and in bow-shocks. In the near future, however, we expect 



that these limitations will be overcome by the instruments 
aboard the Herschel satellite, whose performance in terms 
of sensitivity and of spatial and spectral resolution will be 
closer to those attainable with ground-based observations 
at shorter wavelengths. 

The aims of the present work arc threefold. Firstly, 
we wish to test the ability of a current shock model to 
self-consistently reproduce emission from several differ- 
ent species observed over a wide range of wavelengths. 
Secondly, we will investigate which combinations of these 
observations prove most useful in determining shock con- 
ditions, particularly shock type. Finally, from this analy- 
sis, we intend to identify important diagnostic molecular 
lines, mainly of H 2 0, falling in the wavelength range ob- 
servable with Herschel. As a test case, we have applied our 
analysis to the Herbig-Haro object HH54. This is located 
at the northern edge of the nearby star-forming region 
Chall (d« 200 pc, Hughes & Hartigan, 1992). The object, 
consisting of a complex of several arcsecond-scale bright 
knots, was originally discovered by Schwartz (1977), then 
spectroscopically observed in the optical by Schwartz & 
Dopita (1980) and Graham & Hartigan (1988), who found 
a blue-shifted gas with a radial velocity that becomes in- 
creasingly blue in a long 'streamer' from the north to the 
south and that may be a jet responsible for the excita- 
tion of the HH54 complex. Knee (1992, hereafter K92) 
associates HH54 with a monopolar blue-shifted CO out- 
flow, whose driving source remains unclear after a proper 
motion analysis of the infrared knots (Caratti o Garatti, 
et al. 2006). Molecular hydrogen emission was first re- 
ported by Sandell et al. (1987), while the near-infrared 
spectrum was obtained by Gredel (1994, hereafter G94). 
In the mid-infrared, observations with ISO/CAM (Cabrit, 
et al. 1999), ISO/SWS (Neufeld, Melnick & Harwit, 1998), 
and quite recently, with Spitzer/IRS (Neufeld et al., 2006) 
have shown the H 2 ortho-to-para ratio to lie below the 
equilibrium value expected at the observed gas tempera- 
ture. Far- infrared emission, detected with ISO/LWS has 
been reported by Nisini et al. (1996) and discussed by 
Liseau et al. (1996) in the framework of planar, steady- 
state C-shocks. Conversely, Wilgenbus et al. (2000), who 
compared the predictions of steady-state shock models to 
observations of pure rotational (Neufeld et al., 1998) and 
ro-vibrational (G94) H 2 lines, concluded that no single 
planar C- or J-type shock is able to account for the ob- 
served emission in both pure rotational and ro-vibrational 
lines. In the present work, we will check whether or not 
this result can be confirmed on the basis of a large obser- 
vational data-base. 



The structure of the paper is the following: in Sect. 2 
the observations used for our analysis are presented, while 
in Sect. 3 the adopted shock model is described and ap- 
plied to the data. Our conclusions are then summarized 
in Sect. 4. 
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Fig. 1. H2 2.12/im (continuum-subtracted) Sofl image of 
HH54. The knots of emission are labelled according to the 
nomenclature of G94 and Sandell et al. (1987). The con- 
tours are 3,5,10,15 and 20x the standard deviation to the 
mean background (5xl0 -15 erg s _1 cm -2 arcsec -2 ). The 
orientations of the slit used for spectroscopic observations 
are depicted. 

2. Observations and Results 

Our database consists of data collected with the ESO fa- 
cilities and retrieved from the ISO archive . The relevant 
information on the observational settings is summarized 
as a log in Table 

2.1. Ground-based observations 
2.1.1. Imaging 

A narrow-band image in the H2 1-0S(1) filter was ob- 
tained with Sofl (Lidman, Cuby & Vanzi, 2002) at the 
New Technology Telescope (NTT, La Silla, Chile). The 
observation was conducted by nodding and jittering the 
telescope around the pointed position (a(J2000)=12 ft- 55 m 
50.3 s , <S(J2000) = -76° 56' 23") in the usual ABB'A'mode. 
The raw data were reduced by using standard procedures 
for bad pixel removal, flat fielding and sky subtraction. 
Flux calibration was derived from the observation of a 
photometric standard star during the night. A continuum- 
free image was obtained by subtracting from the narrow 
band image an appropriately scaled K s image, acquired 
just after the observation. The H2, continuum-subtracted 
image is shown in Fig^ where the knots of line emission 
are labelled according to the notation of G94 and Sandell 
et al. (1987). The HH54 complex appears to be composed 
of arcsecond-scale bright knots: those labelled A, B and C 
lie along the axis of the optical "streamer" extending south 

1 available at http://www.iso.vilspa.esa.es/ida/index.html 



Table 3. Medium-resolution spectroscopy parameters. 



Knot 


slit 


a 

V-V cloud 

(km s~ 


Av 


CI 


III 


-15 


19±5 


B 


III 


-12 


13±5 


A 


III 


-12 


16±5 


C2 


IV 


-14 


12±8 


L 


IV 


-13 


20±6 


E 


IV 


-12 


13±7 


K 


IV 


-9 


21±8 



Note: a The estimated uncertainty is within 3 km s 



of HH54; for brevity, we will refer to this axis as the "jet 
axis" in the following. In addition, two further knots (E,K) 
plus (at least) three more diffuse regions are characterized 
by a bow-shaped morphology. Additional knots, roughly 
located along the axis and barely visible in the G94 image, 
have been detected (namely, J,L,M). The (1-0)S(1) line 
flux of the whole HH54 complex is 1.12xl0 -12 erg s -1 
cm~ 2 , i.e. approximately a factor of two and five larger 
than the values measured by G94 and Sandell (1987), re- 
spectively. The brightness variability of the H2 1-0 S(l) 
line on a time scale of about 10 years, pointed out by 
G94, is thus confirmed by the present observations. As 
suggested by Liseau et al. (1996), such variability indi- 
cates that HH54 is a non-equilibrium flow. 

2.1.2. Low-resolution spectroscopy 

Long-slit, near infrared spectra of the HH54 region were 
collected with the ESO facilities during two different runs. 
In 1999, we obtained the 1.55-2.50 /im spectra of the knots 
along the brightest wing of the bow (knots E,K) with 
Sofl, while the jet axis was targetted in 2002 with ISAAC 
(Cuby et al., 2004) at the Very Large Telescope (VLT, 
Paranal, Chile). Four grating settings (SZ,J,SH, SK fil- 
ters) are required with ISAAC to cover the spectrum be- 
tween 0.98-2.50 /im. The adopted slit orientations along 
with the knots encompassed by the slit are depicted in 
Fig'-D The raw spectral images were flat-fielded, sky sub- 
tracted and corrected for the optical distortions along both 
the spatial and spectral directions. Telluric features were 
removed by dividing the extracted spectra by that of a 
blackbody-normalized telluric standard star (spectral type 
09V), once corrected for its intrinsic spectral features. 
Wavelength calibration was derived from the lines of a 
Xenon-Argon lamp, while flux calibration was based on 
the photometry of the H2 1-0S(1) continuum-subtracted 
image, whose associated uncertainty is within 10%. Ratios 
of lines falling in different portions of the ISAAC spec- 
trum are affected by an uncertainty of up to 20%, due to 
the different instrumental responses in the four segments. 
This error was estimated by comparing the fluxes of lines 
present in overlapping parts of the spectrum. Fig.0 and 
Table show the ISAAC spectrum and the line fluxes 
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Table 1. Journal of observations. 



Imaging 



Telescope/Instrument 


Date 


Filter 


FoV 


Pixel Scale 


Exposure Time 








(') x (') 


(7pxi) 


(s) 


NTT/SofI 


6 Jun 1999 


H 2 2.12/zm 


4.9x4.9 


0.29 


300 


ISO/ISOCAM 


20 Oct 1997 


CVF 5.0-16.8A«n 


3x3 


6 


3600 



Spectroscopy 



Telescope /Instrument 


Date 


A 


slit / aperture 


P.A. 


A/AA 


Exposure Time 






(pm) 




(°) 




(s) 


VLT/ISAAC 


11 Jul 2002 


0.98-2.50 


I 


28 


1600 


1600 


NTT/SofI 


6 Jun 1999 


1.53-2.52 


II 


140 


600 


2400 


VLT/ISAAC 


31 Dec 2004 


2.06-2.18 


I/III 


28/160 


8900 


1800 


ISO/LWS 


17 Feb 1996 


43-197 


6>~ 80" 




200 


2740 


ISO/LWS -FP 


23 Sep 1997 


63.13-63.23 


6»~ 80" 




8500 


6000 



detected in knot B, which is among the knots richest in 
lines. Clearly, regions at different degrees of dissociation 
and ionization are encompassed by the slit aperture: low 
excitation conditions are characterized by the presence of 
copious H2 ro-vibrational lines (with v < 4 and a maxi- 
mum excitation energy of ~ 25000 K), while the presence 
of ionized gas at T w 10000 K is evidenced from the detec- 
tion of several forbidden lines from [Fe 11] , [C 1] , [S 11] , [N 1] 
along with hydrogen and helium recombination lines. 



bow-shock (Hartigan, Raymond & Hartmann, 1987). An 
average velocity of « 50 km s _1 along the jet axis, however, 
has been recently estimated from proper motion analysis 
by Caratti o Garatti et al. (2006). We can consider this 
value to be an upper limit to the shock velocity, given the 
possibility that the observed shocks travel in a medium al- 
ready put into motion by previous ejection events: a phe- 
nomenon that is frequently observed in the environments 
of protostellar outflow (e.g. Chrysostomou et al., 2000). 



2.1.3. 2.12 jum medium-resolution spectroscopy 

Medium-resolution spectra of the H2 2.1218 /im line were 
obtained in January 2004 with ISAAC, along the direc- 
tions depicted in Fig.^ (slits I and III). The slit aperture 
is 0.3" x 120", which gives a nominal resolution of about 
8900, i.e. 33.7 km s . The data were acquired and reduced 
with the same techniques outlined in the previous Section. 
Wavelength calibration was performed using bright OH at- 
mospheric lines (Rousselot et al., 2000), which provide an 
accuracy to within 3 km s . The instrumental profile in 
the dispersion direction, as measured from Gaussian fits 
to sky lines, was ~ 31.7 km s^ 1 . We give in TableOHthe 
velocity parameters for all the encompassed knots. The 
peak radial velocity (v-v cloud) was computed with respect 
to the ambient molecular cloud, for which a velocity of 
2.0 kms -1 in the local standard of rest (LSR) has been 
adopted (K92). All the knots appear blue-shifted by about 
10 km s _1 , a circumstance which, although confirming the 
association of HH54 with a blue-shifted outflow, excludes 
the possibility that H2 emission originates in the same gas 
that gives rise to the optical spectrum, where v rac i ranges 
from -20 to -70 km s" 1 (Graham & Hartigan, 1998). The 
observed FWHM is larger than the instrumental width by 
a few km s _1 : the deconvolved values, which give an esti- 
mate of the velocity spread (Aw), range between 13 and 
20 km s _1 . The poor spectral resolution, coupled with the 
low signal-to-noise ratio at the base of the line, prevents 
an accurate measurement of the intrinsic FWZI, which 
could provide some information about the speed of the 



2.2. Space-borne observations 

2.2.1. Mid-infrared spectro-imaging 

Images of HH54 in the H2 pure rotational transitions from 
0-0 S(2) to S(7) were obtained with the infrared camera 
ISOCAM (Cesarsky et al., 1996) aboard ISO, in the guar- 
anteed time CAMFLOWS. Spectral scans were obtained 
with the Circular Variable Filter (CVF) over a wavelength 
range of 16.61 to 5.08 /jm with full spectral sampling, and 
a camera pixel size of 6" yielding a total field of view 
of about 3' by 3'. The raw data files were retrived from 
the ISO archive and reprocessed locally with the CAM 
Interactive Analysis (CIA; Ott et al., 2001) software. The 
data reduction procedure consists of dark current and flat 
field correction and cosmic ray removal. Continuum was 
then removed under each H2 line by fitting a linear base- 
line, and the line flux was obtained by trapezoidal inte- 
gration over 5 spectral elements, yielding individual im- 
ages in each H2 line. A preliminary report of this dataset 
was presented by Cabrit et al. (1999). In the reduced im- 
ages, shown in FigO the same bow-shape morphology de- 
lineated by the 1-0S(1) line is recognizable but the gas 
traced by the rotational lines appears smooth, rather than 
clumpy as in the 2.12 /zm image. This could be ascribed 
to the poor spatial resolution at mid- infrared wavelengths, 
but could also reflect a different mechanism at the origin 
of the excitation of the ro-vibrational and pure rotational 
lines. We will address this point further on in the paper 
(see Sec rO|) . 
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Table 2. Near infrared line fluxes observed in knot B. 



5 



Line id. 


A 


F±AF a 


Line id. 




A 


F± A F a 




M 


(lO^ergcm^s- 1 ) 






H 


(10 _15 crgcm" 2 s" 1 ) 


H2 lines 


2-0 S(9) 


1.053 


0.66±0.07 


3-1 Q(7) 




1.368 


0.4 ±0.1 


2-0 S(7) 


1.064 


0.62±0.04 


1-0 S(9) 




1.688 


0.7 ±0.1 


2-0 S(6) 


1.073 


0.4 ±0.1 


1-0 S(8) d 




1.715 


0.6 ±0.1 


2-0 S(5) 


1.085 


0.69±0.05 


1-0 S(7) e 




1.748 


3.9 ±0.1 


2-0 S(4) 


1.100 


0.3 ±0.2 b 


1-0 S(6) 




1.788 


1.86±0.05 


3-1 S(9)+3-l S(10)+3-lS(ll) 


1.120-1.120-1.121 


1.4 ±0.5 b 


1-0 S 3 




1.958 


15.4 ±0.2 


3-1 S(7) 


1.130 


0.8 ±0.3 b 


2-1 S(4) 




2.004 


0.42±0.09 


2-0 S(2) 


1.138 


0.6 ±0.1 


1-0 S(2) 




2.034 


4.20±0.06 


3-1 S(5) 


1.152 


0.9 ±0.1 


2-1 S(3) 




2.073 


1.63±0.04 


2-0 S(l) 


1.162 


1.2 ±0.1 


1-0 S(l) 




2.122 


9.59±0.05 


3-1 S(3) 


1.186 


0.8 ±0.1 


2-1 S(2) 




2.154 


0.51±0.05 


2-0 S(0)+4-2 S(10) 


1.189-1.190 


0.4 ±0.1 


3-2 S(3) 




2.201 


0.3 ±0.1 


4-2 S(9) 


1.196 


0.4 ±0.1 


1-0 S(0) 




2.223 


2.58±0.05 


4-2 S(8)+4-2 S(ll) 


1.198-1.199 


0.5 ±0.2 b 


2-1 S(l) 




2.248 


1.0 ±0.2 


4-2 S(7) 


1.205 


0.3 ±0.1 


2-1 S(0) 




2.355 


0.32±0.08 


3-1 S(2) 


1.207 


0.4 ±0.1 


3-2 S(l) 




2.386 


0.5 ±0.1 


4-2 S(5) 


1.226 


0.4 ±0.1 


1-0 Q(l) 




2.407 


10.9 ±0.3 


3-1 S(l) 


1.233 


0.5 ±0.1 


1-0 Q(2) 




2.413 


4.5 ±0.4 


2-0 Q(l) 


1.238 


0.55±0.05 


1-0 Q(3) 




2.424 


9.6 ±0.4 


2-0 Q(2)+4-2 S(4) 


1.242-1.242 


0.39±0.04 


1-0 Q(4) 




2.437 


4.2 ±0.5 


2-0 Q(3) c 


1.247 


0.56±0.04 


1-0 Q(5) 




2.455 


3.2 ±0.5 


2-0 0(3) 


1.335 


0.39±0.07 


1-0 Q(6) 




2.476 


2.2 ±0.4 


3-1 Q(5)+4-2S(0) 


1.342-1.342 


0.5±0.2 6 


1-0 Q(7) 




2.500 


6.2 ±0.6 


Ionic lines 


[Ci] 1 ^ - 3 Pi 


0.983 


0.98±0.06 


[Fen] a 4 D 5 / 2 — a 


JJ5/2 


1.295 


1.4 ±0.2 


[Ci] 1 ^ - 3 P 2 


0.985 


2.33±0.05 


[Fen] a 4 D 7/2 — a 


6 D 7/2 


1.321 


2.9 ±0.2 


[Sn] 2 P 3/2 - 2 D 3/2 


1.029 


0.21±0.07 


[Fen] a 4 D 5/2 — a 




1.328 


0.6 ±0.1 


[Sn] 2 P 3 /2 ~ 2 D 5/2 


1.032 


0.30±0.07 


[Fen] a 4 D 5 / 2 — a 


F9/2 


1.534 


1.2 ±0.2 


[Sn] 2 P 1/2 - 2 D 3/2 


1.034 


0.28±0.09 


[Fen] a 4 D 3/2 — a 


F7/2 


1.600 


0.7 ±0.2 


[Sn] 2 P 1/2 - 2 D s/a 


1.037 


0.2 ±0.1 b 


[Fen] a 4 D 7/2 — a 


F9/2 


1.644 


8.2 ±0.1 


Rel2 3 S-2 3 P° 


1.083 


0.53±0.05 


[Fen] a 4 Di/ 2 — a 


4 F 5/2 


1.664 


0.28±0.07 


Pa7 


1.094 


< 0.8 f 


[Fe 11] a 4 D 5 / 2 - a 


%,2 


1.678 


0.58±0.08 


[Fen]a 4 D 7/2 - a 6 D 9/2 


1.257 


8.7 ±0.1 


[Fe 11] a 4 D 3/2 - a 


4 F 3 /2 


1.798 


0.4 ±0.1 


[Feu]a 4 D 3/2 -a 6 D 3/2 


1.279 


1.1 ±0.2 


[Fen] a 4 D 7/2 - a 


4 F 7/2 


1.810 


3.4 ±0.4 


Pa/3 


1.283 


1.5 ±0.4 


Br7 




2.166 


<0.2 / 



Notes: a The uncertainty refers to the rms of the local baseline. 

b signal-to-noise ratio between 2 and 3. 

c blended with [Fen] a 4 D 3/2 - a 6 D 5/2 (A=1.248 /an). 

d blended with [Fen] a 4 D 3/2 - a 4 F 5/2 (A=1.711 /an). 

e blended with [Fell] a 4 P 3/2 - a 4 D 7/2 (A=1.749 /un). 

* 3 sigma upper limit. 



The extent of the regions where emission is detected 
at over 3cr above the sky level, increases from the high- 
est frequency line 0-0S(7) at 5.5 /xm to the 0-0S(2) line at 
12.3/xm, an effect which is due both to the loss of spatial 
resolution at increasing wavelength and to the progres- 
sively lower temperatures (from thousands to hundreds 
of Kelvin) probed by lines with lower excitation energy. 
As one of our aims is to obtain a meaningful comparison 



between far- and near-infrared lines, we have superposed 
the 0-0 contours on the 2.12 /um image, so as to identify 
the regions where the emission overlaps. The ISOCAM 
spatial resolution does not allow us to disentangle struc- 
tures with an angular size of less than 6 arcsec, a dimen- 
sion often larger than that of the near infrared knots. In 
the 5.5 /zm image, which better matches the morphology 
delineated at 2.12/mi, we are able to resolve four peaks 
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Fig. 2. 1.0-2.5 fim low resolution spectrum of HH54B. Lines marked with an asterisk are detected with signal-to-noise 
ratio less than 3. 



roughly corresponding to the three knots from north to 
south (knots C to A) and to the wing (knots E,K). Aiming 
to attribute a mid-infrared flux to each peak, we have 
defined four zones, each centred on the individual peak, 
and containing at least 10% of the peak emission; inside 
each zone the flux has been computed by using the task 
polyphot in the IRAF package 2 . The results are given in 
Table 01 where each zone has been identified by indicat- 
ing the corresponding knots at 2.12/im and the area over 
which the photometry has been computed. Since fluxes 
of mid-infrared lines are poorly affected by the extinction 
(which is in any case low, as will be shown in Sect. 3.1), 
their ratios directly give a precise indication of the lo- 
cal ortho-to-para ratio. Based on ISO/SWS observations, 

2 IRAF (Image Reduction and Analysis Facility) is dis- 
tributed by the National Optical Astronomy Observatories, 
which are operated by AURA, Inc., cooperative agreement 
with the National Science Foundation. 



Neufeld, Melnick & Harwit (1998) estimate a value of 1.2 
averaged over all the region, while Neufeld et al. (2006) 
derive values from 0.4 to 2 in knots C and E,K, which 
they interpret as a legacy of an earlier stage when the 
gas had reached equilibrium at a temperature ^ 90 K. 
Disregarding the ratio involving the 5.5 /^m line, whose 
image seems to be affected by calibrations problems (as 
also suggested by the Boltzmann diagram, see FigE}, we 
measure ortho-to-para ratios of 0.9, 1.3 and 2.3, along the 
jet direction, from knots C to A. These values, which are 
in good agreement with both the estimates of Neufeld, 
Melnick & Harwit (1998) and Neufeld et al. (2006) could 
be interpreted in terms of subsequent shock heating events 
due to an episodic jet. 
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Fig. 3. ISOCAM images of HH54. Contours are in steps of 2a from a 3cr level of 8.3 1CP 16 erg s _1 cm -2 arcsec -1 
(at 5.5 /urn, 6.1 ^m and 6.9/xm), 2.7 1CP 16 erg s -1 cm -2 arcsec -1 (at 8.0 pm), 1.4 1CP 15 erg s -1 cm -2 arcsec -1 (at 
9.6 fim) and 5.5 10~ 16 erg s -1 cm' 2 arcsec -1 (at 12.3 /im). 



Table 4. Photometry of (0-0) lines observed with ISOCAM. 



Term A F±AF (10 -13 erg cm -2 s -1 ) 

C,L B,M A E,K,J 

(/im) (135 arcsec 2 ) (207 arcsec 2 ) (155 arcsec 2 ) (173 arcsec 2 ) 



0-0S(2) 


12.3 


5.5±0.2 


5.5±0.3 


1.4±0.2 


5.9±0.2 


0-0S(3) 


9.6 


3.9±0.3 


6.1±0.4 


3.0±0.3 


7.0±0.4 


0-0S(4) 


8.0 


3.4±0.3 


4.3±0.3 


1.2±0.3 


3.6±0.4 


0-0S(5) 


6.9 


3.6±0.3 


6.5±0.3 


3.1±0.3 


7.7±0.4 


0-0S(6) 


6.1 


1.9±0.3 


3.9±0.4 


1.4±0.3 


3.5±0.4 


0-0S(7) 


5.5 


1.3±0.4 


2.2±0.6 


1.3±0.4 


3.3±0.5 



Note: Due to transient effects at low signal levels, the fluxes could be underestimated by a factor 1.3-1.6 (relative fluxes, 
excitation temperatures, and o/p ratios would remain unaffected). 



2.2.2. Spectroscopy 

The 43-197 /im spectrum of HH54 was obtained with 
the Long Wavelength Spectrometer (LWS) aboard ISO 
and reported by Nisini ct al. (1996) and Liseau et al. 
(1996). The spectrum was taken at the beginning of the 
ISO mission (orbit #92) and originally processed with 
a preliminary version of the Offline Processing Software 
(OLP). Recently, ISO telemetry was reprocessed (with the 
OLP10.1) and the resulting Relative Spectral Response 
Function (RSRF) was higher than in the earlier versions; 
here we report the analysis of this upgraded version of 



the spectrum. The observation was carried out with the 
LWS AOT01 (A/AA - 200), adopting the "fast scanning" 
option and an oversampling of 4 times the spectral res- 
olution clement. The flux calibration, based on Uranus 
observations, has an uncertainty of up 30%, while the ac- 
curacy on the wavelength calibration is about 25% of the 
spectral element (0.29 [im. and 0.60 [im for the short [43-90 
/im] and long [90-197 /im] wavelength range, respectively). 
Post-pipeline processing, which consists of glitch removal, 
averaging of different spectral scans, and correction for 
interference fringes, was performed with the ISO Spectral 



8 T.Giannini et al.: Molecular emission in HH54 

Table 5. Far infrared lines observed in HH54. 



Term 




^vac 


F± A F 




(/xm) 


(pro) 


(10 -13 erg s _1 cm -2 ) 


[Oi] 3 Pi- 3 P 2 


63.18 


63.18 


102± 2 


o-H 2 221-I10 


108.07 6 


108.07 


< 1.2 


CO 23-22 + 0-H2O 4i4-3 03 


113.54 


113.458-113.537 


1.3±0.4 


OH 2 n 3 / 2 ,5/2 — 2 n 3/ / 2 ,3/2 


119.45 


119.23-119.44 


3.0±0.3 


CO 20-19° 


130.40 


130.37 


0.9±0.3 


CO 19-18 


137.20 6 


137.20 


<1 


P-H2O 3l3-2 02 


138.57 


138.53 


1.0±0.3 


CO 18-17 


144.78 c 


144.78 


2.1±0.3 


[0 1] 3 P - 3 Pi 


145.52 c 


145.52 


4.5±0.3 


CO 17-16 


153.19 


153.27 


2.2±0.3 


[Cll] 2 P 3 /2- 2 Pl/2 


157.72 


157.74 


9.0±0.4 


CO 16-15 


162. 73 c 


162.81 


3.9±0.2 


OH 2 n i/ / 2 ,3/2 — 2 ^-1/2,1/2° 


163.04 c 


163.12-163.40 


1.5±0.2 


CO 15-14 


173.63 6 


173.63 


6 ±1 


0-H2O 3 03 -2l2 


174.63'' 


174.63 


2±l d 


0-H2O 2i2-l l 


179.53 


179.53 


7±1 


CO 14-13 


185.85 


186.00 


6±1 



Notes: a Not detected by Nisini et al.(1996); b fixed parameter; c deblended lines; d signal to noise ratio less than 3. 




-100 100 



v - v cioud ( km s_1 ) 

Fig. 4. LWS-FP spectrum of the [Oi]63 [im line. The 

Gaussian fit is shown as a dotted line. 

Analysis Package 3 (1SAP), v. 2.1. In comparison with the 
values reported by Nisini et al. (1996), the new line fluxes, 
given in Table are systematically 30-40% lower. In ad- 
dition, two more lines (CO 20-19 at 130.4 /jm and OH at 
163 /im) have been detected, but we do not observe the 
CO 19-18 line at 137 ^m, which was just detectable in the 
old version of the spectrum (shown by Liseau et al. 1996). 
The low resolution data have been complemented by the 

3 The ISO Spectral Analysis Package ISAP is a 
joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are 
CESR,IAS,IPAC,MPE,RAL and SROAN. 



intermediate resolution spectrum of the [0 1] 63 /jm line 
(see Fig.0J. This was accomplished by coupling the LWS 
with a Fabry-Perot interferometer, thus achieving a reso- 
lution of «8500 (Swinyard et al., 1996). The line, detected 
with a signal-to-noise ratio of about five, peaks at v - v c ioud 
= -12 km s _1 , similar to that of the 2.12 /mi line. Although 
this could suggest a similar origin for oxygen and molec- 
ular hydrogen emission, the 63 /im profile is ss 39 kms -1 
wide, and the corresponding FWZI, computed on the line 
fit by applying the procedure of Hartigan, Raymond & 
Hartmann (1987), is w 80 kms -1 . This value can only 
be considered a rough estimate, given the poor signal-to- 
noise ratio of the detection, nevertheless it indicates that 
the 63 /um emission may arise in a shock with velocity 
larger than (or close to) the H 2 dissociation limit. 

3. Analysis and Discussion 

3.1. Boltzmann diagrams 

Analysis of the H2 line emission gives rise to estimates 
for the reddening towards the encompassed knots and 
the temperature of the emitting gas. These two quanti- 
ties can be simultaneously estimated in the framework of 
the Boltzmann (or excitation) diagram, which is a plot 
of \n(N v j jgj) against E v j/k, where N v j (cm" 2 ) is the 
column density of level (v, J), E v j/k (K) its excitation 
energy and gj = (2 J + 1)(27 + 1) its statistical weight. 
The nuclear spin quantum number is I = 1 for ortho- 
H 2 and I = for para~H 2 . If the gas is thermalized at a 
single temperature, the points in the diagram fall onto a 
straight line: a best fit (by minimization of the deviation 
of the data-points from the line) to the gradient provides 
the reciprocal of the gas temperature. At the same time, 
under the assumption of optically thin emission, the ex- 
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Fig. 5. Boltzmann diagram of the four zones identified in the (0-0) images (see text). In the upper-left panel, sym- 
bols adopted to mark different vibrational levels are shown. In each diagram, the derived extinction, along with the 
temperatures of the vibrational manifolds, are indicated. 



tinction can be evaluated by noting that transitions from 
the same upper level should lie at the same point in the 
diagram (e.g. G94). Thus, having adopted the reddening 
law by Rieke & Lebofsky (1985), we have tuned the Ay 
value to maximize the overlap of data points correspond- 
ing to lines originating in the same upper level. We derive 
Ay = l mag over all of the region. 

In a shock wave there is a range of kinetic tempera- 
tures and the Boltzmann diagram forms a curve. In order, 
therefore, to probe our observations at different tempera- 
tures, columns of the pure rotational lines are plotted in a 
Boltzmann diagram together with those of near-infrared 
ro-vibrational lines. Since the region emitting in 0-0 lines 
has been separated in four distinct zones (as described in 
Section 2.2.1), we have constructed four Boltzmann dia- 
grams, plotted in Fig.0 in each zone, the column densities 
of the ro-vibrational lines have been computed by assum- 



ing an emitting region defined from the 3cx contours of the 
2.12 /im line. The similarity of the four Boltzmann dia- 
grams in Fig. implies that the excitation conditions are 
broadly the same for all of the knots in the HH54 region, at 
least at the spatial resolution (ss ten arcsec) characteris- 
tic of the four zones. Each diagram indicates the presence 
of two temperature components. A 'cold' component at 
Tw 800 K is traced by the 0-0 lines: this is in substan- 
tial agreement with that measured by Neufeld, Melnick & 
Harwit (1998) and Cabrit et al.(1999) on the basis of SWS 
and CAM observations (~ 650 K). All the other lines ap- 
pear thermalized at about ~ 3000 K from north to south 
(knots C to A) and at ~ 2500 K on the wing; noticeably, 
this latter result could simply reflect the lack along the 
wing direction of spectroscopic data in the J band, where 
H2 lines with high excitation energy (i.e. those tracing 
hotter temperature components), are located (Giannini et 
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al., 2002). In the plot, the ortho-to-para ratio was assumed 
equal to the statistical equilibrium value of three: any sig- 
nificant departure from this value should result in a sys- 
tematic displacement of the ortho with respect to the para 
levels. While such behaviour is observed among the v—0 
lines (see also Table , no significant deviations from the 
equilibrium value are appreciable in the data points trac- 
ing the hotter gas components (those with t>>l).This can 
more easily be seen by taking ratios of fluxes of bright 
ortho and para lines belonging to the same manifold (e.g. 
those among 1-0 and 2-1 lines, see Table |2J , for which the 
extinction correction is 20% at worst. 



3.2. Models of Molecular Emission 

In the previous Section we showed that the observed 
emission is quantitatively similar in both the jet and 
the wing of HH54, indicating common excitation mecha- 
nisms and similar physical conditions for the whole region. 
Accordingly, we can examine together the H2 emission 
with the emission of the other molecular species (CO, H2O 
and OH), which are observed with a beam encompass- 
ing the whole HH54 region. Furthermore, we have taken 
the zone encompassing knots B and M (hereafter the "B- 
zone"), in which H2 is emitted over the widest range of 
energies observed, to be representative of the excitation 
conditions in HH54. To explore the diagnostic capabili- 
ties of the observed lines in disentangling different shock 
structures and physical parameters, we have analyzed the 
observations in the context of a full range of shock models, 
including steady-state non-dissociative (C-type) and dis- 
sociative (J-type) shocks and quasi-steady J-type shocks 
with magnetic precursor (denoted C, J and C+J, respec- 
tively, for brevity). For each type of shock we have pro- 
ceeded in the following way: i) the H2 column densities 
in the B-zone were fitted using a shock model, which is 
described in Section 13.2.11 ii) temperature and density 
profiles along with local velocity gradient predicted by 
each best-fitting model were used as input parameters for 
a Large Velocity Gradient (LVG) code to reproduce the 
H2O, CO and OH FIR line emission (see Section l3.2.3|) . 
At this step we can both test the overall consistency of 
the selected models and also identify the FIR lines most 
useful in diagnosing the shock physics; in) finally, in or- 
der to determine if atomic/ionic and molecular emission 
are excited under the same conditions, the O 1, C I, S II and 
Fe II lines fluxes predicted by the best shock model found 
for the molecular lines were compared with the NIR ob- 
servations fSection l3.2.4fl . 



3.2.1. H2 Emission 

The shock code, MHD_VODE, has been described by 
Lc Bourlot et al. (2002) and Flower et al. (2003), while 
the method used to calculate the H2 emission has been 
explained by Giannini et al. (2004) and M c Coey et al. 
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Fig. 7. Velocity and temperature profiles from the models 
of the best-fitting steady-state C-type and J-type shock 
(top and middle panels) and the best-fitting J-type shock 
with a magnetic precursor (bottom panel). The velocity is 
expressed in the shock frame and the pre-shock gas is to 
the left. 



Fig. 8. Top: CO line fluxes as a function of the rotational 
quantum number 3 up . Continuous lines give the fluxes pre- 
dicted by the C-shock model by assuming the collisional 
parameters of McKee (1982) (short-dashed) and of Flower 
(2001) plus Schoier et al. (2005) (long-dashed). The error 
bars represent the la uncertainty. The downward arrow 
indicates a 3cr upper limit; middle: as for the top panel 
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Fig. 9. Top panel: H2O line fluxes as a function of wave- 
length. The LWS spectral range extends up to 200 |xm, but 
we report in the figures all the lines falling in the range 
open to Herschel observations. Data points with the rel- 
ative la uncertainties are represented with open circles, 
while the ortho- and para-H 2 fluxes predicted by the C- 
shock model of Tab. H3 are shown with filled circles and 
triangles, respectively. The downward arrow indicates a 
3ct upper limit. The wavelength range covered by PACS 
and HIFI is indicated, as well. Middle panel: as in the 
top panel but for the steady-state J-shock model. Bottom 
panel: as in the top panel but for the J-shock with a mag- 
netic precursor model. 



(2004). Here, only the aspects most relevant to this study 
are reviewed. 

The code solves for one-dimensional, planar, multi- 
fluid flow and can simulate both steady-state J- and C- 
type shock waves and quasi-steady J-type shocks with 
magnetic precursor. The differential equations that deter- 
mine the abundances of the chemical species are solved 



in parallel with the magnetohydrodynamical conserva- 
tion equations for the neutral, positively and negatively 
charged fluids. 

The calculations for the distribution of population 
among 150 H2 ro-vibrational levels (i.e. up to an energy 
of 3.9xl0 4 K) take into account collisional excitation and 
de-excitation, spontaneous radiative decay, collisional dis- 
sociation and ionization, and reformation of H2 on grains, 
and are solved in parallel with the chemical and dynamical 
conservation equations. This approach is essential to en- 
sure the accuracy of the computed H 2 column densities, 
because the level populations do not respond instanta- 
neously to changes in the physical state of the gas. 

The column density of a level is sensitive to the pre- 
shock density, the shock velocity, the age of the shock and 
the strength of the magnetic field. The response of the 
level populations to these variables, which are parame- 
ters of the code, can be seen in the excitation diagram. If 
the gas is in local thermodynamic equilibrium (LTE), the 
Boltzmann plot is a smooth curve, with little scatter of 
the points about the median line. Departures from LTE 
enhance the scatter. Increases in the shock velocity, v s , 
and the pre-shock density, njj, act to increase the rate of 
collisions that vibrationally excite the H2 molecules and 
thus decrease the departure from LTE. The presence of 
a magnetic field acts to dampen and broaden the shock 
wave; therefore, an increased magnetic field strength en- 
hances the departure from LTE. The magnetic field is as- 
sumed to scale by a factor b with the square root of the gas 
density and is expressed as B(iiG)= &[,rtH(cm~ 3 )] - 5 . The 
assumed dependence of B on nu derives from considera- 
tions of energy equipartition: the magnetic energy density 
is oc B 2 , whereas the thermal energy density is cx n. We 
note that the shock front is assumed to be planar, and the 
parameter b scales the component of the magnetic field 
parallel to the shock front (and perpendicular to the flow 
direction), consonant with the assumed dependence of B 
on nu- Thus, nn = 10 4 cm -3 , 6=1 corresponds to a 
transverse magnetic field strength B = 100 fiG. 
The evolutionary age of the shock wave has a marked ef- 
fect on the excitation diagram. The younger the object, 
the closer is the shock wave to being pure J-type and the 
higher are the column densities of the more highly ex- 
cited ro-vibrational levels. Furthermore, the populations 
of these levels tend to be closer to LTE. In a J-type shock 
with a magnetic precursor, the levels of v — 0, 1 are pop- 
ulated principally in the precursor and exhibit greater de- 
partures from LTE than the higher levels. Furthermore, 
most of the population remains in the v = ground vi- 
brational state due to radiative cascade from excited vi- 
brational states. 

The code was systematically applied over a grid of the 
above parameters for each category of shock (C, J and 
C+J) as follows: 10 3 cm~ 3 < n H < 10 5 cm" 3 ; 1 km s" 1 
< v s < w cr it (in the case of C-type shocks), where w cr it is 
the shock velocity at which a steady-state C-type shock 
becomes discontinuous due to H2 dissociation (Le Bourlot 
ct al., 2002); /iG < B < 500 ^G; and, for ages greater 
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Table 6. Parameters of shock models selected, for each 
category of shocks, as best fits to the H 2 emission. 



Shock type 


^shock 

(km s" 1 ) 


(cm- 3 ) 


Age 

(yr) 


B b 
(A»G) 


C 


37 


10 4 




100 


J 


19 


10 4 




10 


C+J 


18 


10 4 


400 


100 



a n H = n(H) + 2n(H 2 ) 

6 B(jj,G)= &[n H (cm- 3 )] - 5 , b = 1 for the C and C+J shock 
models and 0.1 for the J shock model. 



than that at which a J-type shock with a magnetic precur- 
sor becomes indistinguishable from a J-type shock wave 
(approximately 100 yrs for nn = 10 4 cm -3 and v s = 18 
km s" 1 ). 

From rotational H2 transitions observed towards 
HH54, Neufeld et al. (1998) derived a pre-shock ortho- 
to-para ratio of <0.2. We investigated the effect of the 
initial ortho-to-para ratio on the Boltzmann diagram and 
found that, indeed, a value of <0.2 is required in order to 
reproduce the populations in the v=0 levels. If the max- 
imum neutral temperature is sufficiently high (exceeding 
1300 K for C-type shocks, or 10 4 K forJ-type shocks) then 
the local ortho-to-para ratio attains the statistical equilib- 
rium value of 3 (Wilgenbus et al., 2000). This condition is 
met in the models considered (see Fig. 7) and so the v>0 
levels, which are populated mostly in the hot gas behind 
the shock front, are insensitive to the initial ortho-to-para 
ratio. 

For the steady C and J shock models, we have inte- 
grated line fluxes out to a flow time of 10 5 yr, where the 
post-shock gas has reached thermo-chemical equilibrium. 
For the young C+J shock, the total flow time through the 
shock wave should be less than the time to reach steady 
state for the C precursor, about 10 4 yrs (otherwise one 
would have a pure C-shock, not a C+J). However, the 
exact flow time through the J-cooling zone following the 
C-precursor is not precisely constrained, as it depends on 
the detailed shock history (e.g. Lesaffre et al., 2004a, b), 
and on any time-variability in the inflow conditions. Here 
we have chosen, for illustrative purposes, to truncate our 
flux integration at tj = 3000 yrs, by this time the contri- 
bution to the H2 and ISO FIR lines in negligible. We note 
that that fully time-dependent (and 2D) calculations are 
required for more detailed predictions. 

The first result is that C-shock models are unable to re- 
produce the H2 emission: a J-type component is required 
in order to produce the observed degree of thermalisa- 
tion among the levels. As a demonstration, the excitation 
diagram from the best-fitting steady-state C-type shock 
model (v s = 37 km s -1 , nn = 10 4 cm -3 ) is plotted along 
with the observations in Figure |H1 top panel. This model 
fits the v=l levels reasonably well but overestimates the 



lower energy levels while underestimating the higher en- 
ergy levels. 

On the other hand, two models were found to produce 
a good fit to the H2 emission: a J-type shock with a mag- 
netic precursor of age 400 yrs and y s hock—^ km s" 1 and 
a steady-state J-type shock with v shock — 19 km s . A 
pre-shock density of 10 4 cm -3 is predicted by both mod- 
els. Tableland Fig.[7|report the full set of parameters and 
the velocity and temperature profiles, respectively, while 
the excitation diagrams produced from these models can 
be compared with those derived from observations in Fig. 
|S| The v>l energy levels are fitted well by both models, 
because the J-component, which is largely responsible for 
the population of the higher energy levels in the case of 
the J-type shock with magnetic precursor, is essentially 
the same in each case (see Fig. [S] middle and bottom pan- 
els). The C+J shock model reproduces the population in 
the v=0 levels but underestimates the v=l levels. The 
higher degree of thermalisation arising in the steady-state 
J-type shock results in a greater population in the v=l 
levels but this is at the cost of population in the v=0 lev- 
els and is still insufficient to reproduce the observations, 
a result which agrees with the findings by Wilgenbus et 
al.(2000). 

Therefore, at this step of the analysis, it can be con- 
cluded that a dissociative component is required to fit the 
H2 emission, even if from the H2 data alone, it is not pos- 
sible to distinguish between a steady-state J-type shock 
and a quasi-steady J-type shock with magnetic precursor. 
It is interesting to note, however, that a model of the latter 
type could be used to describe the larger scale geometry 
of the HH54 region. In particular, the smoothness of the 
bow-shape traced by the rotational lines (Fig[3Jl, which 
are excited predominantly in the continuous precursor, in 
comparison with the dumpiness of the ro- vibrational lines, 
which arise in the discontinuity and are more affected by 
the time variability of the driving jet. 

3.2.2. Implications of the derived shock parameters 

We want now briefly to analyze the environmental im- 
pact caused by the physical quantities which derive from 
the shock wave model. Indeed, numerical simulations show 
that the momentum transfer rate from a high Mach num- 
ber jet to the ambient medium has a similar efficiency in 
ID and 3D (Chernin et al., 1994). Therefore, it is mean- 
ingful to compare the momentum transfer rate derived 
from our ID shock parameters with that measured in 
the CO outflow by K92. With our estimated total shock 
area of 60" x 60" (from CO line intensities, see Sec l3.2.3l) . 
a preshock density of njj=10 4 cm -3 and a shock veloc- 
ity of 18 kms -1 , we infer a momentum rate of 3 xlO 
Mq kms -1 yr -1 . In comparison, K92 estimated a lower 
limit of 5 xlO" 5 M Q kms" 1 yr" 1 in the CO outflow, which 
is only a factor of six lower than that derived by our 
shock model. Consequently, our proposed molecular shock 
is compatible with the working surface shock driving the 
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whole CO outflow, although we cannot exclude the pos- 
sibility that other gas components [e.g. un-shocked gas 
(Caratti o Garatti et al., 2006)] may exist and thus con- 
tribute to the momentum transfer balance. 

We can further evaluate the mass-loss rate of the driv- 
ing engine, if the jet speed and the density ratio between 
the jet and the ambient medium (which determines the 
fraction of jet momentum transferred to the cloud; e.g. 
Chernin et al., 1994) are known. The optical "streamer" 
south of HH54, which is presumably tracing the "jet" in 
this system, has a typical radial velocity of -70 km s _1 with 
respect to the cloud (Graham & Hartigan, 1988), hence 
the jet speed is probably at least 100 kms^ 1 . Comparing 
the atomic line ratios in the streamer measured by the 
same authors with the ID models of Hartigan, Morse 
& Raymond (1994) indicates a jet density « 10 3 cm~ 3 . 
Since our H2 shock modelling suggests an ambient molec- 
ular cloud density of 10 4 cm~ 3 , the jet would be strongly 
under-dense, and the efficiency of momentum transfer to 
the HH54 shock would be close to 1. The mass- loss rate in 
the jet would then be < 3 xlO -4 M Q kms -1 yr" 1 / 100 
kins" 1 = 3 xl0~ 6 Moyr- 1 . 

3.2.3. CO, H2O and OH emission 

Aiming to understand the diagnostic capabilities of the 
far-infrared lines in discriminating different shock types, 
we have considered all three models described in Table 
even though steady-state C-type shocks have been already 
disregarded on the basis of the H2 emission fit. Our ID as- 
sumption imposes a major limitation on the prediction of 
line profiles, which can be strongly influenced by the geom- 
etry; hence we will limit our analysis to the computation 
of the brightness of the lines. For each model, the physical 
parameters at each time-step along the shock have been 
iteratively used as input of a LVG code in plane-parallel 
geometry. This gives as output the brightness of the lines 
at each temporal step, which have been summed up to 
obtain the total integrated value to be compared with the 
observations. 

In the LVG approximation, the line ratios are con- 
strained only by the local gas kinetic temperature and 
density, in addition to the ratio of the species number den- 
sity and velocity gradient. Hence, in this part of the anal- 
ysis, we simply compare observed line ratios with those 
determined by the physical conditions provided by the 
shock profiles. The unique free parameters are the pre- 
dicted absolute line brightnesses, whose ratios with the 
observed fluxes constrain for each molecular species the 
angular size of the emitting region, and the temperature 
of the dust, whose continuum emission has the effect of 
pumping the population of the low lying energy levels. 
The latter parameter, however, is relevant only for the 
OH emission fitting (e.g. Thai-Q-Thung et al., 1998, Offer 
& van Dishoeck, 1992). 

The code has been developed for the first 45 rotational 
levels for both ortho- and para-H 2 0, 41 (and 61) levels 



for CO and 24 for OH. Energy levels and radiative de- 
cay rates have been taken from the LAMDA database 
(Leiden Atomic and Molecular Database) 4 . In the H2O 
model, collisions with P-H2 and 0-H2 are considered sep- 
arately in the temperature range 5-140 K (Dubernet & 
Grosjean, 2002; Grojean, Dubernet & Ceccarelli, 2003, 
Phillips, Maluendes & Green, 1996), while at higher tem- 
peratures (up to 2000 K) the collisional rates with He have 
been rescaled to obtain a first order approximation for col- 
lisions with H2 (Green, Maluendes & McLean, 1993). We 
have assumed an ortho-to-para H 2 ratio of 3, as the large 
uncertainties associated with the observations do not al- 
low us to discriminate among different values of this pa- 
rameter. In the case of the OH model we have taken the 
collisional coefficients with 0-H2 and P-H2 given by Offer, 
van Hernert & van Dishoeck (1994) in the range from 15 
to 300 K. Finally, for the CO molecule, we have explored 
the sensitivity of the model to the different collisional rates 
available in the literature. McKee et al. (1982) give the jjq 
coefficients, scaled from the collisional rates with He, to 
compute the downward rates for levels with J up <60 and 
200 K < T < 2000 K: this dataset seems therefore ap- 
propriate for line predictions in high-temperature regime. 
The low-temperature regime, on the contrary, has been 
investigated by Flower (2001) who has computed the col- 
lisional rates in the range 5 to 400 K, including energy 
levels up to 3 up —29 and 3 up =20 for collisions with para- 
H2 and ortho-H2, respectively. Extrapolations to include 
energy levels up to J. up =40 and collisional temperatures 
up to 2000 K have been made by Schoier et al. (2005). 

The results for the CO line predictions are plotted 
in FiglSl the long-dashed line refers to the collisional 
rates from McKee (1982), while the short-dashed one to 
the rates from Flower (2001) and Schoier et al. (2005). 
Major differences (up to one order of magnitude in the 
predicted line brightness) occur for J np 10, i.e. for 
levels mainly populated at the low temperatures. This 
difference is particularly evident in the J-type shock 
model (middle panel) : adopting the coefficients of McKee 
et al. (1982), we obtain a two-component distribution, 
with peaks at J up ~8 and J up ~16, while the coefficients 
of Flower (2001) and Schoier et al. (2005) produce a 
distribution of fluxes with a single peak at J up ~16. Both 
sets of coefficients give similar results in the range of the 
ISO observations, which are rather well fitted by both 
the steady-state C-type shock and by the (C+J) shock. 
On the other hand, the predictions of the steady-state 
J-type model deviate significantly from the observations. 
Therefore, we conclude (having already disregarded the 
C-type shock on the basis of the H2 fit) that the (C+J) 
model is the only one able to account for both H2 and 
CO observations, a result underlying the importance of 
a multi-species analysis in modelling line emission from 
shocks. 



4 available at http://www.strw. leidenuniv.nl~moldata 
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Wavelength" Term E up C J C+J r b 

(pm) (cm -1 ) (ergcm _2 s _1 ) 



75.4 


0- 


3ar 


-2l2 


212.16 


1.1 


10" 


13 


1.1 


10" 


13 


7.2 


10- 


14 


~ 10 


78.7* 


0- 


4 2 3" 


-3l2 


300.36 


5.1 


10" 


14 


3.2 


10" 


14 


1.6 


10" 


14 


<< 1 


82.0* 


0- 


616- 


-5o5 


447.25 


5.4 


10" 


14 


2.9 


10" 


14 


1.3 


10" 


14 


<< 1 


99.5* 


0- 


505' 


-4l4 


325.35 


1.2 


10" 




7.5 


10" 




4.2 


10" 




<< 1 


125.3* 


1> 


4 4- 


-3l3 


222.05 


6.1 


nr 


14 


4.2 


10" 


14 


2.4 


10" 


14 


~ 1 


138.5 


I> 


3l3" 


-2fl2 


142.29 


1.3 


10" 


13 


1.0 


10" 


13 


8.1 


10" 


14 


~ 10 


174.6 


0- 


3o3- 


-2l2 


176.76 


3.2 


10" 


13 


2.4 


10" 


13 


1.5 


10" 


13 


~ 10 


179.5 


0- 


2l2" 


■loi 


79.50 


6.9 


10" 


13 


6.9 


10" 


13 


6.9 


10" 


13 


~ 10 


257.8* 


0- 


3 2 i- 


-3l2 


212.16 


1.9 


10" 


14 


2.5 


10" 


IS 


1.5 


10" 


15 


<< 1 


269.3 


P- 


lor 


-Ooo 


37.14 


2.4 


10" 


13 


3.0 


10" 


13 


2.9 


10" 


13 


~ 5 10 2 


398.6 


I> 


2n- 


-2o2 


95.17 


2.6 


10" 


14 


1.1 


10" 


14 


1.4 


10" 


14 


~ 1 10 2 


538.3* 


0- 


lio- 


-loi 


42.37 


5.6 


10" 


14 


9.6 


10" 


14 


2.0 


10" 


13 


~ 3 10 2 



Notes to the table: a lines whose flux changes by more than a factor of two, depending from the shock type, are marked with 
an asterisk; 6 maximum value of the optical depth along the profile of the (C+J) shock. 



The angular size found from the ratio of the theoret- 
ical line brightness, from the (C+J) model, of the 16-15 
line (at 162 /im) with the observed flux, is about 40", re- 
gardless of the adopted set of coefficients. By considering 
that the LWS aperture (80") represents the Airy disk at 
«100/jm of a point-like source (Swinyard et al., 1996), it 
follows that the derived dimension is a lower limit to the 
true size of the emitting region. This latter, taking into ac- 
count the diffraction effects at 162 /im, can be estimated 
to be approximately 60", i.e. a slightly larger region than 
that traced by the H2 pure rotational emission. The model 
indicates that the far-infrared lines arise mainly from gas 
at T« 500-1000 K, which is in agreement with the tem- 
perature found from consideration of the H2 0-0 lines in 
the Boltzmann diagram. Transitions with 6sS 3 up ^10 (at 
wavelengths between 260 and 433 /im) are predicted to be 
from a few to an order of magnitude brighter than in the 
other two models: their observation would thus provide a 
stringent test to the (C+J) model. This will be possible 
in the near future with the HIFI spectrometer aboard the 
Herschel satellite, but also from the ground with the new 
sub-mm APEX telescope (up to the J=8-7 transition). At 
present, a map of the CO 1-0 emission has been reported 
by K92, from which a flux of w 4 xl0~ 16 ergs -1 cm -2 is 
derived by integrating over « 1 arcmin 2 . Noticeably, this 
value is only a factor of 2.5 lower than that predicted by 
the C+J model, if the coefficients of Flower (2001) and 
Schoier et al. (2005) are adopted. 

Fig.|5|reports the fits to the H2O ortho and para lines. 
All three shock models provide good fits to the observed 
lines, with the marginal exception of the steady-state C- 
type shock model: the poor diagnostic capability of the 
LWS water lines is due to their similar excitation temper- 
atures (of w 50-100 K). To search for lines more sensitive 
to the shock parameters, we have reported in the figure 
the line fluxes predicted for bright H2O lines falling in the 
spectral range covered by the two Herschel spectrometers 



PACS and HIFI. Indeed, under the physical conditions de- 
termined by the considered shock models, some fluxes of 
lines (mainly at the shortest wavelengths) are predicted 
to change by more than a factor of two. These, indicated 
by an asterisk in Tab. will be all observable with a 
signal-to-noise ratio larger than ten with an integration 
time of few minutes with PACS and less than one hour 
with HIFI (Poglitsch, Waelkens & Geis, 2001, de Graauw 
& Helmich, 2001). In the same table, we also report the 
maximum value of the line optical depth along the (C+J) 
shock profile: we note that lines originating in low-lying 
levels (those with E up ^ 200 K) are optically thick, be- 
ing the largest r values found in correspondence of back- 
bone transitions. We also show in Fig. ^| (top panel) the 
brightness and optical depth of the 179 /im line predicted 
by the (C+J) model. The line brightness increases to its 
final value of ~ 10 -4 erg cm -2 s _1 sr" 1 with the sud- 
den increase of the temperature at the J-shock front; a 
corresponding decrease in the optical depth is seen at the 
same point. In the bottom panel, the water column den- 
sity and abundance are plotted: the latter increases at the 
shock front because of the activation of very efficient high- 
temperature reactions (at T^ 300 K), which convert all 
oxygen not locked in CO into water (e.g. Bergin, Neufeld, 
& Melnick, 1998). In the post-shock region, N(H 2 0) lev- 
els out at w 6 10 16 cm" 2 . At the same time, the average 
abundance <x(H20)> over the far-IR emitting region, de- 
termined as the ratio N(H20)/N(H2) integrated over the 
shock up to 3 10 3 yr, is ~ 7 10 -5 . This roughly agrees 
with the estimates by Liseau et al. (1996) and also with 
the values found in a number of outflows from Class 
sources (Giannini, Nisini & Lorenzetti, 2001). Thus, the 
discrepancy between water abundance measurements and 
predictions from classical non-dissociative shock models 
(x(H20) larger than 10~ 4 , Kaufman & Neufeld, 1996) can 
be better reconciled in the framework of J-type shock with 
magnetic precursor models. Moreover, since gaseous wa- 
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13 t 1 1 1 1 — 111 1 11 -8 

100 time(yr) 1000 

Fig. 10. Top: Output of the LVG model (in the frame- 
work of the C+J model): 179 /im line brightness and opti- 
cal depth. The neutral temperature is also shown for com- 
parison. Bottom: Water column density and abundance 
predicted by the (C+J) model. 

ter survives only in high-temperature regions, we expect 
that its emission arises in compact parts of the shock: by 
normalizing the predicted 179 /jm line brightness to the 
observed flux, we derive, even considering the diffraction 
at the line wavelength, 9 «30 arcsec. Therefore, the water 
emission region is more compact than that emitting the 
warm CO by about a factor of four. Noticeably, a differ- 
ence in the shock area for CO and H2O lines, signals some 
inconsistencies in ID geometry, since in this approxima- 
tion the gradients within the cooling zone (perpendicular 
to the shock surface) are already taken into account and 
thus the emitting area should be exactly the same for all 
the species. 

As regards the OH molecule, we have detected only 
two lines at 119 and 163 /im, each connecting the lowest 
energy levels of the two ladders ( 2 n 3 / 2 and 2 n 1 / 2 ) into 
which the OH energy level diagram is split. These lines 
have a similar excitation temperature and critical density, 
thus not allowing us from discriminating between differ- 
ent shock categories. Using the C+J shock structure, the 
LVG computation predicts a 119/163 ratio of 14, which is 
substantially different from the observed value of 2. This 
discrepancy could be due to the pumping effect of a dust 
continuum, the heating of which is not accounted for in 
our shock model. Indeed, the 163 /im line brightness is in- 
creased more than that of the 119 /im by dust pumping 
(Offer & van Dishoeck, 1992). 

In summary, our multi-species analysis leads us to con- 
sider the (C+J) shock as the model which best reproduces 



Table 8. Molecular parameters derived for the J-type 
shock with magnetic precursor model. 



Species 


N 


9 


<x> 




(cm" 2 ) 


(arcsec) 




H 2 


7.4 10 20 






CO 


1.2 10 17 


60 


1.6 10~ 4 


H 2 


5.4 10 1B 


30 


7.3 10" 5 


OH 


1.3 10 14 




1.8 10" 7 



Table 9. [O i]63,145/im brightness as predicted by shock 
models. 



Model 


1(63) 


1(145) 




(erg cm 2 


s" 1 sr" 1 ) 


C 


9.9 10 -7 


9.1 10" 8 


J 


3.2 10~ 8 


1.5 10^ 9 


C+J 


1.3 10 -6 


1.4 10" 7 



the overall observed molecular emission. The parameters 
relative to the FIR emission are given in Tab. |H1 Here, we 
report the molecular column densities and the fractional 
abundances predicted by the (C+J) shock model and in- 
tegrated along the shock profile, along with the angular 
sizes of the emitting regions derived by the ratio of the 
predicted brightnesses with the observed fluxes. 

3.2.4. Atomic emission 

Of the atomic and ionic species detected in HH54, oxygen 
is the only one involved in the same chemical reactions 
as H2O and OH. Therefore, it seems reasonable to expect 
that a sizeable fraction of the emission at 63 and 145 /tm 
should arise in the same shock as the molecular lines. In 
Table we give the line brightness of the [O i] 63 and 
145 /tm predicted by the three shock models considered in 
this paper. Even for the highest values (predicted by the 
(C+J) model), an emitting area fifty times larger than the 
ISO beam is required in order to reproduce the observed 
flux of the 63 /im line. This suggests that oxygen emission 
originates in gas excited by a different mechanism to that 
which gives rise to the H 2 line emission. Both the FWHM 
of the 63 /im line profile and the observed 63/145 ratio 
of 23 are typical of radiative J-shock environments (e.g. 
Hollenbach & McKee, 1989), which could form at the apex 
of the bow-shock or in reverse shocks along the jet. If we 
compute the mass-flux from the 63 /tm flux using the for- 
mula Mshock = M /yrx (L(63 /im)/L Q ) (e.g. Liseau 
et al., 1997) we find M shock = 1.2 1CT 6 M Q /yr, which is 
similar to the jet mass flux estimated in Sec l3.2.2l 

With regard to the other atomic and ionic species, none 
of the observed lines can be accounted for by the (C+J) 
shock model: the predicted intensities of [Fen], [Ci] and 
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[S n] are approximately 4 orders of magnitude smaller than 
are observed. Indeed, Spitzer maps (Neufeld et al. 2006) 
have shown that the bulk of the fine structure emission 
(from [Neil], [Si], [Sin], [Fen]) comes from a defined re- 
gion close to the head of the bow: our results suggest that 
the degree of dissociation and ionization in this region are 
much higher than in the gas in which the H2 emission orig- 
inates. The same authors state that the fine structure lines 
originate in fast J-type shocks of velocity « 35-90 kms -1 , 
in agreement with the (rough) estimate we derived from 
the [0 1] 63 /jm line for the shock velocity. 
In conclusion, the atomic/ionic emission appears to trace 
a different (faster and more ionizing) shock than the H2 
emission. Similar results have already obtained by us when 
fitting the near- infrared lines observed in HH99 (M c Coey 
et al., 2004): in that case, in a self-consistent scenario, 
two J-shock components, one with a magnetic precursor 
and one fully dissociative, were required to account for 
the H2 and the ionic emissions, respectively. However, we 
postpone a full analysis of this topic to a later paper. 



4. Conclusions 

We have presented a multi-frequency /multi-species anal- 
ysis of the molecular emission in HH54, by means of one of 
the most complete spectroscopic data-base in the infrared 
(from 1 to 200 /xm) ever collected for an HH object. The 
H 2 lines (1-12 (im) coming from levels v=0 to v=A with 
excitation energies between 2000 and 25000 K have been 
interpreted in the context of a state-of-art shock code. Two 
models (a steady-state J-type shock and a quasi-steady J- 
typc shock with magnetic precursor) yield satisfactory fits 
to the data, allowing us to discard steady-state C-type 
shocks as a mechanism for the line excitation. The out- 
put parameters (temperature and density profiles, molec- 
ular species number density and velocity gradient) from 
these models are adopted as input to an LVG computa- 
tion to interpret the emission of CO, H2O and OH in the 
far-infrared. This second step of the analysis allows us to 
definitively select the J-type shock with magnetic precur- 
sor as the one able both to account for the overall emis- 
sion of H2, CO and H2O. None of the considered models, 
however, are able to reproduce the observed OH emission, 
which probably could be accounted for by considering the 
heating of dust grains in the shock model. A further, fully 
dissociative, high-velocity shock is thought to be respon- 
sible for the observed atomic/ionic lines. If it is the case, 
a picture may emerge in which the optical and near-IR 
ionic lines trace a fragmented jet shock, the near-IR H 2 
lines trace ambient J-shocks driven by the jet shocks, and 
the v=0 H 2 lines trace the continuos precursor propagat- 
ing into the ambient cloud. 

Finally, the adopted approach has been extended to pro- 
vide predictions of the best diagnostic H 2 line intensities 
in the range accessible to the forthcoming facilities HIFI 
and PACS aboard the Herschcl satellite. 
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